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ABSTRACT 


As devices shrink, the current density through interconnects increases proportionally 
making new materials a necessity for industry growth. Carbon nanofiber (CNF) and 
carbon nanotube’s (CNT) potential for high current density make them a possible 
replacement for metal contacts. Learning the limitations of CNFs and CNTs is important 
if they are to be used in next-generation electronics. As current density increases, heat is 
generated throughout the CNF structure. This heating eventually leads to breakdown as 
the temperature reaches the bonding energy of the Carbon-Carbon (C-C) bond, the bond 
between two carbon atoms. The resultant reaction is the vaporization of the carbon, 
eliminating electromigration problems common with metal interconnects. 

The physics of breakdown of CNFs is poorly understood. The CNF interconnects’ 
heating under a voltage sweep between two electrodes is modeled in this thesis. A 
working model was created with Silvaco ATLAS using experimental data provided by 
Santa Clara University (SCU). An analytical solution was found for the heat generation 
occurring within the device. The simulation does not show the breakdown occurring; 
instead, it accurately predicts the temperature and electrical characteristics of the device. 
This model will aid in the analysis of CNFs’ reliability and potential future integration 
into the next generation electronics. 
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EXECUTIVE SUMMARY 


The purpose of this work was to develop a device simulation that would accurately model 
the maximum current density and the heating within a carbon nanofiber (CNF) 
interconnect. A CNF is a type of carbon nanotube (CNT) that is one possible material 
that may replace copper interconnects in the near future. As such, it is necessary to 
develop a method for analyzing their capabilities. A successful model would have the 
ability to identify potential shortcomings and advantages of CNFs as compared to current 
interconnect technology. Experimentation done by a group at Santa Clara University 
(SCU) was the basis for developing the technology computer aided design (TCAD) 
simulation. 

The experiments from the SCU group analyzed various breakdown and transport 
characteristics of CNFs. Unfortunately, they have no method for monitoring the 
conditions of the device during operation. This group has performed mathematical 
modeling on their structures; however, the functionality of those models is limited. No 
simulation has previously been developed to accurately model a CNF device under 
operation. The Silvaco TCAD products are utilized in this thesis to create a model for 
prediction of current density and heating effects inside the device. Developing a model 
that can accurately simulate a CNF allows for faster testing methods and the ability to 
identify the root causes of failure mechanics. 

The device was developed to test CNFs of various lengths under an applied 
voltage sweep. The device characteristics the SCU group developed during their 
experimental testing were used to define the thermal and contact model parameters. 

Using knowledge of CNTs, we incorporated other physics-based characteristics about the 
device into the model. 

CNFs and CNTs are unlike other interconnects because of their geometry. The 
cylindrical structure of CNFs makes them more difficult to model. As a result, the best 
approximation of the three-dimensional (3D) structure was converted into a two- 
dimensional (2D) cross section of the active regions of the device. Multiple length 
interconnect structures were created for testing. Two main sets of simulations were 


XV 



performed. One set examined the dependenee of maximum eurrent density on the length 
of intereonneet. The other examined the effeet of the CNF coming into contact with the 
substrate material at the center point of the interconnect length. A third simulation 
compares the heating effects in the 2D cross-section model versus a 3D model of the 
same device. 

Simulations of the structures showed that the model could closely approximate 
the data achieved experimentally. For longer length interconnects, the data matches the 
slope of the fitted data from SCU experimentation. Possible causes for the discrepancies 
are either poorly defined physics of the material or a result of the 2D geometry or both. 

While this model may have limited functionality, it provides a starting point for 
more complex modeling of CNF materials. As the model matures, it would be possible 
to apply any electrical parameter to the contact to measure the materials response. 
Utilization of this model could validate the functionality and reliability of CNFs for use 
in interconnects. 
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I. 


INTRODUCTION 


Interconnects are the simplest and the most common components in electronic 
circuits and, as a result, there are a wide variety of types. According to the International 
Technology Roadmap for Semiconductors (ITRS), interconnects are one of a list of 
factors that need to be considered for the future of electronics. The ITRS differentiates 
interconnects into two categories based on gate length of transistors. These categories are 
split between the near term and high power applications, greater than 16 nm, and the long 
term and low power applications, less than 16 nm. As the connection distance decreases 
different factors become important due to various quantum effects. Interconnects that fall 
in the near term/high power applications category are considered in this thesis. 
According to the ITRS, we see that there are five critical factors that contribute to 
interconnect engineering in this group. They include the type of material used, reliability, 
metrology, manufacturing integration, and manufacturing efficiency. The first three 
factors are addressed in this project. 

Electronic interconnects currently utilize copper, but there is a strong push to 
replace them with materials with greater conductivity capabilities. Carbon nano tubes 
(CNTs) offer these properties as well as the possibility for greater reliability. 
Researchers at Santa Clara University (SCU) work with CNTs with slightly different 
properties called carbon nano fibers (CNFs). CNFs are CNTs with non-tubelike structures 
within their outer shells. SCU has been studying the thermal and transport characteristics 
of CNFs for the purpose of integrating them into the next generation devices. Because 
the ability to manufacture these devices has current limitations, the researchers at SCU 
are only able to make one test sample at a time, which leads to a shortage of usable 
samples for extensive testing. In actuality, a large number of CNFs are produced when 
they are grown [1]. The problem is that applying them to a device requires precision, and 
there is currently no automated process to accomplish this task. 

The goal of this project was to recreate some of the experiments done by SCU on 
these CNFs used as interconnects and create a working computer model that could 
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predict and analyze the structure. The model was used to analyze maximum current 
density with respect to structural changes to the CNF interconnect. The results were then 
used to predict the failure limits of the CNF interconnect material. The use of physics 
based computer modeling for the study of new material for use as interconnects is 
demonstrated in this thesis. 

A. RESEARCH APPROACH 

A semiconductor simulation program to construct a physical model of a CNF 
interconnect and simulate its operation under various conditions to determine thermal 
characteristics during operation as well as maximum current density was used in this 
thesis. The program utilized was Silvaco, which is a suite of programs that work together 
to perform physics-based simulation studies of semiconductor devices. The primary 
programs used for this project include DevEdit, DevEditSD, ATEAS, TonyPlotSD, and 
TonyPlot. There were multiple iterations to the program in order to create a working 
model. 

The first phase of the project involved creating two-dimensional (2D) and three- 
dimensional (3D) models of the interconnect device using DevEdit and DevEditSD, 
respectively. DevEdit is a program that can graphically define a device, and DevEditSD 
is the same except that a third axis of a structure can be defined. DevEdit was the only 
software tool that allowed the construction of a circular region, which was necessary for 
the 3D CNE structure. The only way to get the proper dimensions for the 3D device was 
to draw a cross-sectional view of the device using DevEdit3D. By defining the sizes of 
the regions in the third dimension, the structure was generated. The mesh, defining 
where the simulator will solve equations, was also created within DevEdit. Once this 
iteration was complete, the structure could be viewed in TonyPlot and TonyPlot3D. 

In the second phase of the project, the correct materials for each region were 
defined so the structure would have characteristics similar to its real-world counterpart. 
All of the material characteristics for the CNE section were modified from values 
originally associated with silicon carbide (SiC). These changes were all applied in 
ATEAS. 
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The third phase involved verifying that the structure and new material parameters 
could produce results similar to the SCU group. This involved applying a voltage sweep 
and measuring the current density and temperature. In this final phase, some of the 
factors examined included length of the CNF and the effect of supporting a portion of the 
CNF with the substrate. Although the length of the CNF and contacts were changed, the 
substrate of the device was kept constant to maintain proper scale. 

The completed simulations were exported to TonyPlotSD and TonyPlot to analyze 
the heating in specific sections of the CNF. Since there were simulations from multiple 
lengths of the CNF, the simulation had to be run individually for each length. The 
resultant data was used to investigate the dependence of current density and thermal 
characteristics versus various structural factors of the device. 

B. RESEARCH OBJECTIVE 

The goal of this thesis was to develop a computer model of a CNF interconnect to 
investigate the thermal characteristics of a CNF material under various electrical stress. 
There have been papers published that deal with length and contact material for a DC 
mode [2-9]. With a computer simulation, it is possible to introduce a variety of input 
characteristics while maintaining a constant structure. Structural inconsistencies between 
nanotubes can be a problem that limits testing capability. Because of the size of the 
devices, it is incredibly difficult to get the kind of detailed measurements that can be 
generated from a computer model. The research done in [2-11] was only able to 
determine the breakdown current of the CNFs with respect to length and some basic 
thermal characteristics. The 2D simulation strived to create something that could closely 
predict these breakdown characteristics. 

C. IMPORTANCE 

As with any new technology, it is important to establish limits of operation. 
Copper is becoming difficult to manage as device dimensions shrink. This means the 
future of interconnects depends upon new materials and novel ideas. CNTs are one of the 
most promising to replace copper. Although CNFs are not currently used in 
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manufacturing, there will be a push for intereonneets with higher current density 
eapabilities. Beeause they are at least one order of magnitude larger than CNTs, CNFs 
ean be a good replaeement for larger scale deviees that require high current density. 

D. THESIS LAYOUT 

This thesis is divided into five ehapters. Baekground relevant to working with 
CNTs and intereonneets is provided in the next ehapter. The SCU experiments and the 
development of the models are detailed in Chapter III. In Chapter IV, the experiments 
are detailed and analyzed. Finally, in Chapter V, the eonelusions drawn from the 
experiments are explained and reeommendations for improvements and future projeets 
are provided. 
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II. BACKGROUND 


The three most important things to understand for this project are the many types 
of carbon nanotube-like structures (allotropes), interconnect technology, and previous 
simulations of similar structures. To best understand how this element works, it is 
necessary to examine the simplest structures and work towards the more complex. After 
analyzing the properties of the various structures, one can see how this material can be 
incorporated into current interconnect technology. 

A. CARBON 

Electrical and thermal conductance, mechanical strength, and weight-to-strength 
ratio are just some of the reasons as to why carbon and, more specifically, CNTs and 
graphene are of such great interest to the scientific community. CNT’s official discovery 
occurred in 1991 by Sumio lijima, a Nippon Electric Company (NEC) researcher at the 
time [12]. The first CNTs were needle-like tubes composed of multiple layers. lijima 
had created the first experimentally proven, multi-walled carbon nano tube (MWCNT) 
[2]. He later discovered single-walled carbon nanotubes (SWCNT) in 1993 [13]. Since 
lijima’s discovery, study of CNTs has increased enormously. This has led to more 
research being conducted on the manufacturability of CNTs and the effects of replacing a 
multitude of current electronic materials with CNTs. 

Carbon is the only element that has isomers from zero dimensions (OD) to three 
dimensions [1]. Even though all the structures are composed of the same element, their 
properties can vary simply based on their structure. The simplest and perhaps easiest 
isomer to create is graphene. CNTs are rolled up graphene sheets. CNTs come in two 
different categories, SWCNTs and MWCNTs. Easily, you get the CNE which is a 
MWCNT with variable structures on the interior of the innermost tube. This is the order 
in which the different structures will be covered. 
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1. Graphene Properties 


a. Structure 


Graphite is the most common allotrope of carbon composed of multiple 
layers of graphene. Graphene consists of a repeating pattern of hexagonal shapes on a 
planar surface, where each point of the hexagon corresponds to a carbon atom. This 
means that any carbon atom in the structure has three neighboring carbon atoms to which 
it is bonded. In Figure 1, the unit cell and Brillouin zone are shown along with unit 
vectors [1], 


a\ 
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2 2 
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where a is defined as the lattice constant of the carbon structure. 


( 2 ) 


In Figure Ib, the triangle FMK, shown by the dotted line, determines the region 
being solved for the energy dispersion relations. These points are important for 
determining physical characteristics later. 


Unit Cell Brillouin zone 



Figure 1. (a) The unit cell and (b) Brillouin zone are defined for a graphene sheet. 

After [1] 
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b. Physics 

2 

A carbon atom contains four bonds, three a-bonds that hybridize in a sp 
configuration and a 2pz orbital [1]. The a-bond is the strongest form of covalent bond. 
These kinds of bonds are formed by direct overlapping between atomic orbitals. The 2pz 
orbital is perpendicular to the graphene plane, creating a n covalent bond. A 71 -bond is an 
indirect covalent bond. These bonds usually occur between two barbell-shaped parallel p 
orbital shells. The bonding is indirect because the shared electron clouds do not overlap 
the nucleus of the atom and, as a result, are much weaker [1]. 

The three a bonds are used to bond the atom to adjacent carbon atoms. As a 
result, the n bond becomes the most important for determining the solid state properties 
[1]. The energy dispersion relations for the 71 -bond and the 71 * anti-bonding band are 
illustrated in Figure 2. The points in Figure 2 correspond to the Brillouin points seen in 
Figure lb. 



Figure 2. The energy dispersion relations for 2D graphite are shown for the 
Brillouin zone. The inset graphs the energy dispersion along the high symmetry 
directions of the triangle TMK shown in Figure 1. After [1]. 


According to Saito et al [1], since there are two n electrons per unit cell, these two 
71 electrons fully occupy the lower n band. A detailed calculation shows that the density 
of states at the Fermi level is zero. This means that two-dimensional graphite is a zero- 
gap semiconductor. The existence of a zero gap at the K points gives rise to quantum 
effects in the electronic structure of carbon nano tubes [1]. A simple approximation for the 
electronic structure of a graphene layer is provided by [1] as: 
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The variables kx and ky eorrespond to the x and y direetions of the waveveetor k. 
Equation 3 is used to obtain a simple approximation for the eleetronie dispersion relations 
for earbon nano tubes. This dispersion relation is very important for ealeulating the 
energy dispersion and eleetronie strueture for CNTs. 

2. CNT Properties 
a. Structure 

As its name suggests, earbon nanotubes are tube-like struetures eomposed 
of earbon atoms or, more speeifieally, graphene sheets of earbon atoms. CNTs eome in 
two different varieties, SWCNT and MWCNT. The SWCNT has only one layer of 
graphene rolled into a tube-like shape, whereas MWCNT have multiple layers. The 
number of layers ereates a drastie differenee in the properties and physies of the 
structures [1]. 

CNTs are elassified into three eategories based on their symmetry: armehair, 
zigzag, and ehiral. The first two elassifieations are aehiral, whieh means that any 
nanotube with these struetures has a mirror image with an identieal strueture to the 
original tube. The method for elassifying nanotubes involves determining its ehiral 
veetor Ch- The parameter Ch can be expressed by the unit veetors uj and a 2 of the 
hexagonal lattiee defined as [1]: 

C;, = na^ + ma 2 = {n, m), (n, m are integers, 0 < |m| < n). (4) 

The absolute value of Ch, L, eorresponds to the eireumferenee of a eross seetion of 
a nanotube. The variables n and m are used to define what kind of tube has been made. 
When n = m, an armehair nanotube is ereated. A zigzag nanotube is ereated when m = 0. 
Any other eombination of n and m will generate a ehiral nanotube. The differenee 
between these struetures ean be seen in Figure 3. 
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Figure 3. Classification of carbon nanotubes: (a) armchair, (b) zigzag, and (c) chiral 

nanotubes. From [1], 


Once Ch is defined, the diameter of the nano tube di can be calculated from [1] 

d^=LI n, L=|C^| = -sJCi^ -C^ =ay]rd +nd +nm. (5) 

Another important vector that defines a nanotube is the translational vector T. 
The vector T is normal to the chiral vector and is defined to be the unit vector of a one¬ 
dimensional (ID) carbon nano tube by [1]: 

T = tjaj +t 2 a 2 = (where integers). (6) 


The variables ti and C are defined as [1]: 


2m+ n 2n + m 

l/t- n Ct- D 


(V) 


where dR is the greatest common divisor of {2m + n) and {2n + m). With the vectors T 
and Ch, the ID unit cell can be defined. When the vector product operator is applied to 
the two vectors and divided by the area of a hexagon, the number of hexagons per unit 
cell N is obtained defined by [1]: 

IC^xTI 2{m^ +n^ +nm) 21} 

A = P^ = - ’- = ^. (8) 

l^j X ^ 2 1 d,^ a di^ 

Thus, there are 2N carbon atoms in each unit cell of the carbon nanotube [1]. 
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Figure 4. The unrolled honeycomb lattice of a nanotube with Ch = (4, 2). From [14], 

How vectors T and Ch would look on a graphene sheet is illustrated in Figure 4. 
When points O and A are connected and B and B' are connected, a nanotube can be 
constructed. The resultant nanotube is defined as Ch = (4, 2). This is determined by 
using the unit vectors a\ and ai. 

Though there are other parameters that are used to define carbon nanotubes, most 
of them are not important for basic understanding. Table 1 is a complete list of all the 
parameters of CNTs. However, once Ch and T are determined, all the other values can be 
calculated from them. 

All of these parameters play a role when determining the physics of an individual 
CNT. The parameters that are the most important are the n and m values and the 
resulting diameter of the CNT. With these two parameters, the electrical characteristics 
of the CNT can be calculated. 
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Table 1. 


Parameters for carbon nano tubes are shown. From [14]. 


Symbol 

Name 

Formula 

Value 

°c-c 

carbon-carbon distance 


1.421 A (graphite) 

a 

length of unit vector 

VSflcc 

2.46 A 

fli. a? 

unit vectors 


in (Jr.y) coordinates 

b,.b. 

reciprocal lattice vectors 


in (Jt.y) coordinates 

c* 

L 

chiral vector 

circumfeicnce of nanotube 

Cff nt, + msj * 
i = |C*| = aV/t^ nm 

/f,m: integers 

0 s |m| s n 

d, 

diameter of nanolube 

i. 4- + nm 

d. — — ^ - Q 

w r 


6 

chiral angle 

yim 

sinP = — - 

2v/i^ + + nm 

0^ |e| < 30° 



2n + m 

cos 0 = ,- - - 

2yn'' + nm 

•iim 

Ian 0 = - - 

2rt + m 


J 

the highest common divisor of (/t.rn) 



dK 

the highest common divisor of 
(2rt + /«, 2m + n) 

^ (d if /I — m not a multiple of 3d 

" (3d if n - m a multiple of 3d. 


T 

translational vector of 1D unit cell 

1 — /[»( + /j*) = 

2m + n 

Ok 

ti.tj: integers 



2/f + m 

Or 


T 

length of T 

It 


N 

iiuinlier of hexagons pei 1D unit cell 

2(n’ + + nm) 

w 

Ok 

2A' • «( /unit cell 

R 

symmetry vcclort 

R = pa, + q%i = ( p,q) 
d = mp — nq. 0 s p s n/d, 0 s i? < m/d 

p^q". integerst 

M 

number of 2v revolutions 

M ~ ((2n 4- n})p 4 (2/n + f\)Q\/dft 

NH. « 2 V/C* + (iJ 

A/: integer 

R 

basic symmetry operation i 

R = (lf|r) 



rotation operation 

, , -W / 4‘L\ 

2x/ 

t! : radians 

7 

translation operation 

II 

T.x: length 


^{p.q) arc uniquely determined by d = mp - nq, ■jubjert to conditions staled in table, except for ^Kxag tubes for which 
= (n.O). and wc define p = \. q = —I, which gives M = 1. 

}K and R refer to the same symmetry operation. 
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The physics behind carbon nanotubes can get quite complicated and, in 
some cases, are far more involved than is needed for this project. The physics that are 
necessary for effectively defining this project are covered in this section. There are two 
main points that must be covered. The first important point is the method for determining 
if a CNT is metallic or semi-conductive. The second involves defining the general 
electrical structure of CNTs. 

One of the unique properties of CNTs is that they can be either metallic or semi- 
conductive depending on their structure. The effect of chirality in determining if the 
nanotube is metallic or semi-conductive is illustrated in Figure 5. 



Figure 5. The electronic structure of a CNT given the chirality vector. 

One third of the vectors possible produce metallic nanotubes. 

The remaining nanotube combinations are semi-conductive. From [14]. 

This ratio of metallic to semi-conductive nanotubes becomes especially important 
when dealing with MWCNT. Since about one third of the cylinders of a MWCNT are 
conducting, certain electronic properties, such as the electrical conductivity of nanotubes, 
will be dominated by the contributions from the conducting constituents, and the non¬ 
conducting constituents will play almost no role [1]. This is a very crucial point for this 





Because this project deals exclusively with MWCNT, only the electronic structure 
of metallic nanotubes is covered here. As such, it is important to know how and why a 
metallic CNT forms. The condition for metallic energy bands in a CNT is illustrated in 
Figure 6a. 

(b) 

K 

M 

K 

Figure 6. (a) The condition for metallic energy bands is shown, (b) The multiple 

Brillioun zones of a CNT are shown. From [1]. 

The vectors Kj and K 2 are defined as [1]: 

A; = ^ (-C^i + 6^2 )> ^2 = ^ ("^^1 - ”^2 )> (9) 

where b\ and 62 are defined in equation 2 and t\ and C are defined in equation 7. The 
term N is defined in equation 8 and also defines the number of Brillouin zones. In Figure 
6 a, the line WW' corresponds to the first Brillouin zone of this one-dimensional material. 
The vector K 2 defines the length of the WW' line. The Ki vector defines the distance 
between subsequent Brillouin zones up to N iterations. The many Brillouin zones of a 
carbon nanotube with a C/, = (4, 2) is illustrated in Figure 6 b. 

The Brillouin zone is important for determining whether or not the nanotube has 
metallic properties. According to the energy dispersion of the 71 -bands shown in Figure 2, 
the Brillouin zone has zero band gap at the vertexes of the hexagon. In Figure 6a, those 
same points correspond to K and fC. 
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The condition for obtaining a metallic energy band is that the ratio of the length of 
the vector YK to that of Ki in Figure 6a is an integer. Since the vector YK is given by 
[ 1 ]: 




( 10 ) 


the condition for metallic nano tubes is that {2n + m) or equivalently i^n-m) is a multiple 


of three [1]. This is the reason why one third of nanotubes end up being metallic. Once a 
nanotube is determined to be metallic or semi-conductive, the electronic structure of the 
CNT can be defined. 


The electronic structure of metallic CNTs is fairly straightforward. At finite 
temperatures, the band gap of the CNT is zero. This means that very little energy is 
required to excite electrons into the conduction band. Also, for all metallic nanotubes, 
independent of their chirality and diameter, it follows that the density of states per unit 
length along the nano tube axis is a constant given by [1]; 




8 


( 11 ) 


where a is the lattice constant of the graphene layer and |t| is the nearest neighbor Carbon- 


Carbon (C-C) tight binding overlap energy [1]. This means there are a finite number of 
states for conduction. 


An important characteristic pertaining to semiconducting nanotubes is the energy 
gap. The energy gap of semiconducting nantubes depends upon the reciprocal nano tube 
diameter d, [1], 




ha. 


c-c 


d. 


( 12 ) 


independent of the chiral angle of the semiconducting nanotube, where c = a / >/3 is 
the nearest-neighbor C-C distance on a graphene sheet [1]. This means that as the 
diameter of the nanotube increases, the band gap of the CNT decreases. Therefore, since 
this project deals with MWCNT and large diameter CNTs, the band gap of devices is 
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very small or zero even for the semi-conductive nano tube shells. All of these properties 
are the same for CNFs with some slight differences. 

3. CNF Properties 

CNFs are MWNTs that are slightly more complicated than CNTs in that they 
have large diameters and variable internal structures. CNFs are between 25-200 nm, 
which is one to two orders of magnitude greater than most SWCNT and MWCNT. From 
equation 12, these large diameters imply that the bandgap for the semi-conductive CNT 
shells will be at or near zero. Any nano tube with a diameter greater than 14 nm has a 
bandgap comparable to the thermal activation energy [1]. 

Aside from the multiple outer walls, CNFs have internal structures composed of 
failed nano tubes. Failed nano tubes are graphene shells that might be at an inverted angle, 
so that instead of creating a tube, a cone is generated. In some cases, a tube-like structure 
is created, but it does span the entire length of the nanotube. This kind of structure is 
refered to as bamboo-like. The difference between a regular nanotube and the internal 
nanofiber structure is illustrated in Figure 7. When the term a=0, the structure maintains 
a tube-like structure, but as a increases, cone shaped structures begin to form. 



Figure 7. Comparison of nanotube and nanofiber structures is shown. From [15]. 
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In normal MWCNT, the interior of the structure is empty, which means there is 
no conduction through this region. CNFs do have these interior structures but, because 
the distance between the layers is so large, electrons still primarily conduct along the 
exterior of the CNF layers. There can be thermal conductivity through this interior 
region that will contribute to the overall thermal coefficient. The thermal features and 
pertinent physics of the CNF structures are covered in more detail when the SCU 
experiments are discussed. Another important feature of CNFs is that they are easier to 
create then the more ordered CNTs. 

4. CNT/CNF Synthesis Methods 

As mentioned earlier, the first nanotubes were discovered by accident by Sumio 
lijima as a byproduct of one of his experiments [12]. Carbon nano tubes can be grown in 
a variety of ways. The method for growth determines what kind of CNTs are created. 
Some methods are more effective at producing particular kinds of CNTs. Since they 
were first discovered, many new methods have been developed. Some of these methods 
include chemical vapor deposition, arc discharge, and laser ablation. 

a. Chemical Vapor Deposition (CVD) 

CVD is the most widely studied method because of its ability to produce 
CNTs continuously. CVD is the method used by the SCU group to make their CNFs. If 
the proper conditions are defined, this method could be a very good way to synthesize 
large quantities of CNTs under controlled conditions [1]. 

The growth process for CNTs involves the use of a catalyst metal to begin the 
formation of CNT structures. Some of the metals used are nickel (Ni), iron (Fe), and 
cobalt (Co). In order to grow CNTs, any one of these metals are patterned on a substrate 
and placed in a reaction chamber. Then the chamber is heated up to about 1100 C, and 
hydrocarbon gases are pumped in. The gases react with the metal catalysts, slowly 
depositing carbon atoms with a hexagonal lattice. The leftover hydrogen atoms are 
flushed from the chamber by more incoming hydrocarbons. The catalization process and 
the synthesis of CNTs are illustrated in Figure 8. 
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The resulting nanofibers are primarily MWCNT with some SWCNT present. Due 
to their low energy formation, the struetures have poor erystallinity until they are heat- 
treated. There have been improvements to this method that involve the use of plasma to 
facilitate catalyzation. 



Catalyst: Ni, Co, Fe . 

Figure 8. CVD catalyzation process is shown. From [16]. 
b. Arc Discharge 

The arc discharge method uses a different approach to generate CNTs. 
Typical conditions for operating a carbon arc for the synthesis of CNTs include the use of 
carbon rod electrodes separated by a small distance with a voltage across the electrodes 
and a dc electric current flowing between the electrodes [1]. Carbon begins to deposit on 
the negative electrode. This method has the ability to make both MWCNTs and 
SWCNTs. No catalyst is needed to create the MWCNT. SWCNTs can be grown using 
similar catalyst to those used for the CVD method. This method creates limited 
quantities based on the carbon electrode used since the carbon electrode transforms into 
the nanotubes. 


c. Laser Vaporization 

Laser vaporization is an efficient method for synthesizing bundles of 
SWCNTs. This method is performed by placing a Co-Ni/graphite target in a heated flow 
tube, and then two sequenced laser pulses are used to evaporate the target. Flowing 
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argon gas removes the entrained nanotubes from the high temperature region and deposits 
them into a water-eooled Cu eolleetor downstream. By varying growth parameters, the 
average nanotube diameter and diameter distribution ean be varied [1]. 

B. INTERCONNECT TECHNOLOGY 

1. Current Interconnects 

Currently, Cu is the primary material used for intereonneets in deviees. Cu has 
been in use for over a deeade as intereonneets, but new problems are emerging as the 
industry moves forward. During a period in whieh the reliability of Cu-based 
intereonneets is expeeted to rapidly deerease, the reliability requirements per unit length 
is rapidly inereasing [16]. As deviees shrink, new materials are needed to overeome the 
problems assoeiated with eurrent metal intereonneets. The problems most eommonly 
eneountered are maximum eurrent density, eleetromigration, resistanee, and signal delay. 

As deviees shrink, the eurrent density through intereonneets inereases 
proportionally, making new materials a neeessity for industry growth. The maximum 
eurrent densities reeorded for Cu intereonneets are on the order of 10^”10^ AJcvcv. The 
eopper intereonneet is now routinely used with minimum feature sizes down to 130 nm 
[17]. When intereonneet dimensions reaeh the nanometer range, the maximum eurrent 
densities do not allow for large potential differenees between eontaets and breakdown 
oeeurs mueh more quiekly. 

Another faetor that ean eause intereonneets to fail is eleetromigration. 
Eleetromigration is the transport of atoms from one physieal loeation to another on metal 
intereonneets. As atoms move from one point on an intereonneet to another point, 
thinning oeeurs at one end of the deviee, and eurrent density in that region inereases 
rapidly. As intereonneet dimensions deerease, the eurrent density assoeiated with a given 
voltage inereases, eausing an eight-fold inerease in the eleetromigration rate with every 
halving of the line width (for a fixed thiekness-to-width ratio) [16]. 

Resistivity in Cu intereonneets is also a problem as feature dimensions shrink. 
The seattering of eleetrons by surfaees and grain boundaries leads to inereased resistivity 
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as line widths are decreased [16]. This resistance becomes much more noticeable for 
long interconnects. As a result of the resistance increase, there is much greater signal 
delay in large scale devices. The proposed alternative to Cu interconnects is CNT and 
CNF interconnects. 

2. CNF Advantages 

CNTs and CNFs have numerous advantages over metal interconnects, making 
them an important material to consider as a replacement. In nearly all the above 
categories, CNTs outperform Cu interconnects. 

One of the most important selling points of CNTs is that they have the potential 
for much higher current densities (10*~10'° AJcvci) [17-19]. These densities are two to 
three orders of magnitude greater than those found in Cu. A contributing reason for this 
difference is that a CNT’s mean free path is in the micrometer range compared to the ~40 
nm mean free path of Cu [20]. The large mean free path allows ballistic transport 
through the CNT. Ballistic transport results in reduced resistivity and strong atomic 
bonds, which provide tolerance to electromigration [20]. 

Electromigration occurs differently in CNT interconnects than it does in metal 
interconnects. Instead of metal from the interconnect material moving from one site to 
another, a CNT interconnect transports metal between the contact points of the CNT. 
Electromigration occuring on a CNT is illustrated in Eigure 9. Even though 
electromigration still occurs, it is greatly reduced as a result of the ballistic transport of 
the CNT. 



Eigure 9. Electromigration of metal on a CNT is shown. 
Atoms move along the surface of the CNT. Erom [16]. 
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In addition, the resistivity of the CNT is smaller beeause of the ballistie transport. 
This leads to long CNTs offering signifieantly redueed resistanee eompared to Cu [16]. 
For the most part, in shorter length intereonneets the resistanees are eomparable. 

Beeause of inereased potential eurrent density and deereased eleetromigration 
effeets and resistanee, CNTs ean be used in smaller deviee struetures with fewer negative 
effeets. Overall, CNTs and CNFs provide a possible alternative to Cu for future 
intereonneet teehnology. 

C. PREVIOUS WORK 

There have been numerous groups that have made models for CNTs. Only one is 
eovered in this work. The model eovered provides an equivalent eireuit of a metallie 
MWCNT intereonneet. With this model eertain assumptions ean be made about the basie 
eleetrieal eharaeteristies of a metallie MWCNT and, therefore, CNFs. 



Figure 10. Simplified equivalent eireuit of a metallie MWCNT intereonneet. From 

[20]. 

The work by Srivastava et al [20] derives an equivalent eireuit for metallie 
MWCNTs. In their model they eonsider faetors sueh as the number of walls and eontaet 
resistanees. The simplified eireuit is illustrated in Figure 10 [20]. 

The terms Ce and Cq eorrespond to series and parallel eapaeitanees, respeetively, 
between shells. The term Rc is the eontaet resistanee, and the term M is the number of 
shells. The kinetie induetanee per unit length Lk, quantum eapaeitanee per unit length 
Cq, and resistanee R are the terms that define the eleetrieal properties of any given shell. 
This model simplifies further into Figure 11 [20]. 
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Figure 11. Simple equivalent circuit model of a metallic MWCNT interconnect. From 

[20]. 


When M is large, Cq can be neglected, which means the capacitances of a 
MWCNT interconnect is smaller than that of a SWCNT. In addition, the resistance and 
inductance of all metallic shells are parallel and M times smaller than that of a SWCNT 
[20]. Using this information, we see that a CNF also has small values for its resistance, 
capacitance, and inductance [20]. 


D. CHAPTER SUMMARY 

In this chapter, carbon materials and structures, interconnects, and previous work 
modeling CNTs were examined. It was determined that carbon has many structural 
parameters that determine if the CNT/CNF is metallic or semi-conductive. Some of these 
parameters are dependent on the growth methods used. It was found that interconnect 
technology could greatly benefit from the addition of CNTs as a replacement for current 
Cu interconnects. The previous work examined another group who has created models of 
CNT interconnects. 

In the next chapter, the SCU experiments are analyzed. Also, an introduction to 
the Silvaco simulation software and its implementation for creating the test structures are 
examined. 
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III. DEVICE AND SIMULATION 


In this chapter, the proeedure for ereating the simulation is eovered. The first step 
to aeeomplishing this was to define the parameters for the experiment. As diseussed 
previously, the work done by the group at SCU was the inspiration for ereating the model 
and was used extensively to define the CNF parameters of the deviee. In the following 
seetion, the SCU experiments are detailed and the information extraeted was defined. 

The seeond part of this ehapter provides details about the simulation software and 
the generated model. The basies of Silvaeo’s programs are eovered along with a 
deseription of the proeess for ereating struetures and defining deviee properties. Then, 
using data from the SCU experiments, test parameters were input and the model was 
finalized for testing and verifieation. 

A. SCU EXPERIMENTS 

1. Overview 

The primary method of verifieation for this projeet was based on the work done 
by Santa Clara University [2-11]. They eondueted numerous experiments on basie 
earbon nano fiber intereonneets to determine their breakdown eharaeteristies. Some of the 
eharaeteristies that they studied inelude temperature dependenee, the eomposition of the 
eontaets, length dependenee, transport phenomena, and the effeet on the CNF while in 
eontaet with the substrate. Multiple struetures were implemented in order to test all of 
these parameters. 

The basie deviee is eomposed of a silieon dioxide (Si02) substrate with two gold 
eontaets separated by a distanee of mierons. The two gold eontaets are about 100 nm 
thiek and eonneeted by a CNF. The lengths of the CNFs vary between 1.5-8.6 pm and 
have diameters of between 100-200 nm [2]. The CNFs in these experiments are grown 
using the plasma-enhaneed ehemieal vapor deposition method with Ni eatalyst. All of 
the deviees were measured to their breakdown limit by eurrent stressing. The primary 
test parameter studied was the eurrent density through the deviee beeause it is elosely 
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related to the amount of joule heating being generated. The four structures that were 
tested by the group at SCU are illustrated in Figure 12. The Au contact by itself with 
suspended and supported structures, respectively, are illustrated in Figure 12a and 12b. 
The W-Au contact with suspended and supported structures, respectively, are illustrated 
in Figure 12c and 12d. For each set, a scanning electron microscope (SEM) image and a 
simple structural representation are displayed. 


(a) Suspended (b) Supported 



CNF CNF 



Current suppty Cuaent supply 


(c) Suspended with W deposition (d) Supported with W depositcn 



Current supply Current supply 


Figure 12. Test structures and SEM images from SCEl experiments. Erom [2]. 
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In some of their experiments, Yamada et al. [2-6] examined the effeets of laying 
the CNF so that it eame in eontaet with the substrate, as shown in Figure 12b and 12d. 
The CNF ean either be eompletely suspended over the SiOi substrate or partially in 
eontaet with it. The term S is defined as the ratio of the length in eontaet with the Si02 
substrate to the CNF intereleetrode length. When S = 0, the CNF is eompletely 
suspended, and when S = 1, it is eompletely supported by the substrate. Beeause the 
primary method of failure is assumed to be eaused by Joule heating, heat ean diffuse into 
the substrate when the CNF eomes in eontaet with the substrate. This eontaet inereases 
the resulting eurrent density possible through the CNF strueture. When breakdown 
happens, it was proposed that oxidation of the earbon oeeurs [5]. 

Another faetor that the group at SCU experimented with was the effeet of 
different eontaets on the eontaet resistanee and thermal breakdown of the CNF 
intereonneets. Contaet resistanee was a major faetor in a number of the papers. The 
reason for this is that Au has diffieulty bonding to the CNF, ereating large resistanee 
values. In order to deerease this eontaet resistanee, the group used eurrent annealing. 
Current annealing is the proeess of applying a set amount of eurrent through the 
intereonneet, with the intention of eliminating defeets. The initial resistanee at the 
eontaet is on the order of MQ. It was found that after about three rounds of the eurrent 
annealing proeess, the eontaet resistanee bottomed out with a resistanee in the kQ range 
[7]. The W-Au eontaet intereonneets did not benefit from this eurrent annealing; 
however, the addition of the W ereated a better initial eontaet, with a resistanee in the kD 
range. For this projeet, only the W-Au eontaet intereonneets were examined. The next 
seetion eovers the test data that was used to ereate and verify the output of the simulation. 

2. Test Data 

The simulations in this thesis aehieved an approximate mateh to the test data 
shown in Figure 13. The maximum eurrent density versus inverse intereonneet length is 
illustrated in Figure 13. There are three different measurement sets based on the strueture 
at the eontaets: suspended CNFs with Au eleetrodes (solid eireles), and suspended (solid 
triangles) or Si02-supported (open triangles) CNFs with W-Au eleetrodes [2]. 
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(13) 


The fitted eurves are based on the one-dimensional heat transport equation [3-5], 


d^AT{x) 2 AT’ ut 2 u 1 

-r- a AT = -bJ , b = —, 

dx K(J 

where J is the current density and a is the dissipation factor measuring the effectiveness 
of the heat dissipation to the contact material with units of inverse length. The termziris 
the local CNF temperature at x measured from the ambient temperature. Finally, the term 
b equals l/(Ka), where k is the thermal conductivity and a is the electrical conductivity 
[3]. In equation 13, the left half of the equation, heat diffusion and heat dissipation, 
generally depends on location, whereas the right side, heat generation, does not. The 
temperature of the devices wre estimated through a thermogravimetric measurement [5]. 




Figure 13. Measured maximum current density as a function of reciprocal CNF length 

is shown. From [3, 5]. 


The next section incorporates the experimental data from these experiments into 
an ATLAS simulation. 
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B. CNF ATLAS SIMULATION 


I. Silvaco Simulation Software 

The simulation software suite used for this projeet was provided by Silvaeo 
International. The software suite is eomprised of many different programs that are 
designed for simulating eleetronie deviees. The programs used for this projeet from the 
software suite inelude ATLAS, Deekbuild, DevEdit, DevEditSD, TonyPlot, and 
TonyPlotSD. 

The ATEAS deviee simulator was used for the modeling of the CNE intereonneet. 
ATEAS is a speeialized programming language that allows a user to ereate a physieal 
model of an eleetronie deviee. The user ean define the struetural parameters and physies 
of the model. ATEAS allows up to 20,000 nodes and 100,000 nodes when solving 2D 
and 3D simulations, respeetively. A node is a point on the strueture of the deviee where 
the differential equations for eleetrieal and thermal models are solved and eonvergenee is 
tested for the solutions. All of these nodes together ereate a mesh. The mesh for this 
work was defined using the DevEdit program and is deseribed later. 

Convergenee of solutions ean be a eonstant problem when performing 
simulations. Some methods to reduee the possibility of problems are to optimize the 
meshing, determine the best step size when biasing, and to use the proper numerie 
method. ATEAS has the eapability of using multiple methods for obtaining solutions. 
The Newton solution method obtains solutions for eoupled equations, whereas the 
Gummel method solves the deeoupled equations. The Bloek method is a hybrid of the 
two. The Bloek method solves faster than Newton, but may not obtain a solution that is 
as aeeurate as Newton method [21]. The user has the ability to ehoose a eombination of 
the methods to inerease the ehanees of eonvergenee. Eor this projeet, the Newton and 
Bloek methods were used for the 2D simulations, and a hybrid solver ealled Halfimplieit 
was used for the 3D simulation. 

When solving for partieular parameters in deviees, ATEAS has groups of 
equations that are aetivated. The group of equations that are used for temperature 
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calculations are defined as GIGA. ATLAS GIGA self-healing simulator incorporates 
models for mobility, self-heating effects, and lattice temperature [21], ATLAS has built- 
in default values for the material parameters. The user has the ability to define all of the 
parameters present in the equations to accurately simulate the material’s thermal 
properties. 

Deckbuild is a text editor and runtime environment that allows the user to run and 
modify ATLAS programs. The other programs used are described as they are 
encountered. 

2. Device Construction 

In order to accurately construct the device, the programs DevEdit and DevEditSD 
were utilized. Both programs provide the ability to draw cross-sections of the device. In 
DevEdit, the user has the ability to define polygon shapes and map them out on a 
coordinate plane. DevEdit also allows the user to define the mesh and can automatically 
generate a mesh for the structure. 



Eigure 14. Basic DevEdit model (2.0 pm). 
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Figure 15. DevEdit Region menu sidebar. 

In Figure 14, the basie window for both DevEdit programs ean be seen with one 
of the deviee struetures drawn. It should be noted that the origin of the graph is at the top 
left eorner, whieh makes the deviee seem like it is upside down with the aetive region at 
the bottom of the figure. At the time of this projeet, both DevEdit and DevEdit3D only 
funetion on a maehine running Red Hat Einux OS. 

The eommands for making a new strueture are very straightforward. The first 
step is to ereate a new region by using the region drop down menu. Onee seleeted, a new 
sidebar pops up that allows the user to give a name to the region, define the material, and 
ereate the shape of the region. Figure 15 is provided as a reference for defining a region. 

The shape of the region can be created two different ways. The first method is 
accomplished by inputting the points of a polygon or shape underneath the polygon box 
on the right sidebar. The user needs to remember to hit the enter key, so that the program 
registers the change to any point values. For the second method, the user simply draws 
the shape of the region by clicking where on the coordinate plane they wish to place the 
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points. The program automatically generates lines between points and defines the 
structure in the polygon window in the right sidebar. 

Once all the regions are in place, the mesh for the structure can be created. The 
user has the option of using the built in meshing ability of DevEdit or manually defining 
the mesh structure after the structure is brought into ATLAS. The built-in meshing 
ability was used for this project. Special attention was given to the auto mesh constructed 
for the device because it was necessary for making the device converge and accurately 
simulating the device physics. 

The particular structures developed are described next, including all the region 
definitions and the mesh construction. 

a. 2D Structure 

For the purpose of this project, two distinct drawings were created. The 
2D models needed to be drawn as a cross-section of the CNF splitting it along its length. 
This device is shown in Figure 14. The large section in the figure is the substrate 
material. If this were a real device, the substrate would be two orders of magnitude 
thicker. This difference in substrate thickness does impact the thermal properties and is 
addressed later under the thermal contact condition part of this chapter. 

The other materials in Figure 14 are the gold and tungsten contacts and the CNF 
materials. The important thing to note about the contacts is that the tungsten is required 
to test this device in 2D. Without it, the multiple materials that the CNF is composed of 
are not connected, and the outer layer would contribute nothing to the simulation. Since 
the tungsten contact is required for the 2D simulations, only structures with W-Au 
contacts were tested. 

The CNF structure was defined by three separate regions. There were two regions 
of material that made up the outer walls of the CNF. Between these two regions was an 
insulating layer of material that corresponds to the interior structure of the CNF. An 
insulating layer was chosen because there was no electrical conduction through the center 
of the CNF, but there was thermal conduction taking place. The lengths of these regions 
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were always equivalent for every length of CNF tested. Since there is no material that 
currently has the properties of a CNF, 4H-SiC was chosen as the material to modify for 
the exterior and Oxynitride was chosen for the interior. The important parameters of the 
material were changed, so it did not matter which materials were used. 

Multiple structures were generated with the primary focus being on creating 
longer interconnect lengths. This involved shrinking the gold contacts while maintaining 
the same amount of contact with the CNF structure. As a result, the length of the CNF 
materials increase as the gap between the gold contacts increases. A total of six different 
2D structures were generated for the length dependence tests. The lengths were chosen 
so as to best fit the testing data from SCU. 

Another testing group was created by using three of the above structures and 
adding in a block of substrate material that connected the center of the CNF to the main 
substrate. The TonyPlot software is used to illustrate a structure with 25% of the span 
supported by the substrate in Figure 16. The three different lengths were tested with 25% 
and 50% of their length supported. The contact was always placed at the center point of 
the CNF. A complete set of structures generated is shown in Table 2. 

DEVEDIT 

Data from 2D500L25S.str 
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Figure 16. TonyPlot image of five micron interconnect meshed with 25% contact 

with substrate. 
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Table 2. List of 2D structures created for simulations. 


Length (pm) 

1.25 

1.75 

2.00 

2.50 

3.75 

5.00 

Suspended 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Supported (25%) 

Yes 

No 

No 

Yes 

No 

Yes 

Supported (50%) 

Yes 

No 

No 

Yes 

No 

Yes 


The mesh for these structures was generated by the built in meshing function. In 
order to refine the mesh in the places of greatest interest, new mesh constraints needed to 
be defined. Mesh constraint definitions determine the distance between solve points. 
The user can specify constraints for particular groups of materials or individual materials. 
Only the maximum height and maximum width were modified to obtain mesh definition 
for this project. The mesh was focused on the CNF materials and contacts. All the mesh 
constraints applied to the 2D structures are contained in Table 3. A smaller number in 
Table 3 corresponds to a more defined region of the structure. 

Table 3. Mesh constraints for the 2D structures are shown. 


Material 

Max. Height (pm) 

Max. Width (pm) 

Si 02 (substrate) 

0.15 

0.15 

Gold (Au) 

0.025 

0.05 

Tungsten (W) 

0.025 

0.05 

4H-SiC (CNF) 

0.005 

0.005 

Oxynitride (CNF) 

0.005 

0.005 


When importing the DevEdit file into ATLAS, the width of the device needed to 
be determined. Even though the structure is technically 2D, in order for the simulation to 
solve properly, a third dimension needs to be defined for the device. Eor these 
simulations, the default width of one micrometer was chosen. It should be noted that the 
resulting structure looks nothing like a CNT/CNE. The CNE layers will look like 
rectangular blocks. However, since the primary characteristic being measured was 
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current density, the geometry of the simulated device was not as important as the material 
characteristics. For the most accurate structural representation of the device, a 3D model 
was created. 


b. 3D Structure 

A cross-sectional picture of the 3D structure needed to be drawn to create 
the radial dimension of the CNF. The cross-section for the 3D device is illustrated in 
Figure 17. It should be noted that DevEdit is the only program in the Silvaco suite of 
products that allows the creation of circular regions in conjunction with rectangular 
regions. 

Any circle or arc created using DevEdit is actually a polygon with a set number of 
sides. The program allows the user to define how accurately they want the curved region 
to resemble a circle. Eor this experiment, the greatest accuracy was used, which places a 
point at every degree for a total of 360 points. The drawback to choosing this highly 
defined structure was that the mesh structure needed to be extremely refined around the 
boundary points of the CNF structure. This circular pattern was used for both regions of 
the CNF section. The internal region was composed of an insulator with a diameter of 50 
nm, and the outer layer had a thickness of 25 nm, thus giving the CNF a diameter of 100 
nm. 



Figure 17. Cross-section of 3D device with circular CNF. 
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The key difference between DevEdit and DevEditSD is the ability to define the 
third dimension of a device. This is done by defining the z-dimension to various depths. 
One thing to note is that if the user closes the program and then goes back to edit the 
structure, all of the z- dimension values must be input again. This issue also applies to 
the mesh constraints defined for either DevEdit programs. 

The 3D model was not as extensively tested as the 2D simulations. There was 
only one model created for the purpose of determining the heating characteristics in a 
geometrically accurate CNE. 

3. Material Properties 

The primary focus of this project was to design a simulation that could 
approximate the qualities of a CNE material. In order to accomplish this, a known 
material needed to be modified to have the same characteristics found in a CNE. The 
parameters that were changed included the band, recombination, and thermal parameters. 
All changes to the characteristics occur within the ATEAS code. 

The material chosen to be modified was 4H-SiC. It was chosen because 4H-SiC 
is not as common as other materials such as Si within the ATEAS code. As a result of 
extensive testing, some materials, such as Si, defined in ATEAS have a more refined set 
of characteristics that if modified can generate errors in the simulation. Since insulators 
have very few parameters that contribute to a simulation. Oxynitride was chosen as the 
interior insulating material. 

The band parameters of the 4H-SiC material were the greatest contributing factor 

to the material characteristics. As discussed previously, CNEs are MWCNTs with an 

irregular internal structure. Since MWCNT are always metallic, the CNE materials had 

to be metallic. However, as a result of the imperfect formation process, electron traps 

occur in the structure causing the CNE to exhibit semi-conductive qualities. In order to 

account for this property, the two conducting regions were defined as semiconductors 

with a band gap equal to the thermal activation energy (~26 mV). With such a small 

band gap, the material still acted like a conductor but allowed the user to have greater 

control over the physics in the material. This control was very important for measuring 
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the current density. If the material had been defined as a conductor, the ATLAS 
simulation would not have produced current density as an output. 

The other band parameters that played a role were the density of states and the 
permittivity. The density of states has enormous effects on the breakdown of the device. 
As stated earlier in the text, metallic CNTs have a finite density of states. Therefore, the 
density of states for the CNF materials was set at a constant value according to equation 
11. The permittivity of both CNF materials was estimated based on the permittivity of 
similar carbon materials. For this simulation, the permittivity of diamond was used. 

Although there may be some recombination occurring due to traps, the primary 
heating method for the device was from Joule heating. Instead of trying to account for an 
unknown number of traps, the recombination parameters in the device were all set to 
zero. This might have caused some slight error in the finished simulation; however, any 
error should have been negligible compared to the Joule heating effect. 

The thermal parameters were taken directly from the experimental data. ATLAS 
calculates heat based on the lattice heat flow equation [21] 

dT 

C^ = V(xVrj + Lr, (14) 

dt 

where C is the heat capacity per unit volume, k is the thermal conductivity, H is the heat 
generation, and is the total lattice temperature. In general, the thermal conductivity is 
temperature dependant. ATLAS acknowledges this by having multiple models that are 
available for use. According to the SCU data, the thermal and heat conductivity of the 
CNF material are constant. The data was unclear about which portion of the CNF the 
conductivity data applied to. As a result, both the active and passive regions of the CNF 
were defined with these conductivity values. 
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Table 4. 


Passive region (Oxynitride) material parameters ehanged are shown. 


Type of Parameter 

ATLAS identifier 

Changed Value 

Band Parameter 

Permittivity 

5.4 

Thermal Parameter 

te.eonst 

0.12 


HC.A 

1.75 


HC.B, HC.C, HC.D 

0 


Table 5. Aetive region (4H-SiC) material parameters ehanged are shown. 


Type of Parameter 

ATLAS identifier 

Changed Value 

Band Parameter 

Permittivity 

5.4 


EG300 

0.026 


NC300, NV300 

3x10'^ 


NC.F, NV.F 

0 

Reeombination 

taunO, taupO 

0 

Parameters 

nsrhn, nsrhp 

0 


ksrhtn, ksrhtp 

0 


ksrhen, ksrhep 

0 


ksrhgn, ksrhgp 

0 


augn, augp 

0 


augkn, augkp 

0 


kaugen, kaugep 

0 


kaugdn, kaugdp 

0 


Etrap 

0 


Copt 

0 

Thermal Parameters 

te.eonst 

0.12 


HC.A 

1.75 


HC.B, HC.C, HC.D 

0 
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A complete list of the parameters and their ehanged values are shown in Tables 4 
and 5. The band parameters affeet the band strueture of the deviee. The term 
permittivity defines the permittivity of the deviee. The term EG300 defines the bandgap 
of the strueture. The rest of the terms in the band parameters seetion of Table 5 eontrol 
the density of states. The initial density of states in the eonduetion and valenee band are 
defined by the terms NC300 and NV300, respeetively. The terms NC.F and NV.F are an 
exponential faetor that modifies the density of states with respeet to temperature. These 
were ehosen to be zero beeause the density of states needed to remain eonstant. 

The thermal parameters HC.A, HC.B, HC.C, and HC.D affeet the heat eonduetion 
of the deviee and te.eonst speeifies a eonstant thermal eonduetivity value. 

The reeombination parameters in Table 5 are all zero and are not defined beeause they do 
not eontribute to the simulation. If traps were being eonsidered as another possible heat 
souree, the terms etrap, taunO, and taupO would be given values. The general eode used 
to implement these ehanges into the model is eontained in the Appendix. 

4. Thermal Contact Conditions 

The Si02 layer on the bottom of the deviee was mueh thieker than the CNF and 
eontaet layers. The substrate of the deviee was two orders of magnitude thieker than the 
deviee modeled. Sinee this exeess material only makes a thermal eontribution to this 
model, it ean be replaeed by an equivalent thermal eontaet. 

Within the ATFAS program, it is possible to define a thermal eontaet in the 
thermoeontaet statement with the parameter a. From the ATFAS User’s Manual [21], a 
is determined from the thermal resistanee of the material as shown by, 

« = (15) 

a 2S 

where k is the thermal eonduetivity and S is the surfaee length from a point on the 
interfaee to some radius. The term a is measured in W/em^-K. Equation 15 implies that 
the thermal resistanee is determined by the ratio of the width to the eross-seetional area. 
Sinee the thermal eonduetivity of Si02 is known, only the geometry of the represented 
region was needed. 
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The geometry of the substrate region was calculated by using the three- 
dimensional model, where the substrate surface was 20 gm and the depth of the substrate 
was 200 gm. The thermal resistance of a substrate with this geometry was calculated to 
be 1400 W/cm -K. The contact was placed at the bottom 100 nm of the substrate. It was 
assumed that in all cases the same amount of excess substrate would be present. 
Therefore, this value was used for both the 2D and 3D simulations. The implementation 
of this contact can be seen in the Appendix. 

5. Simulation Parameters 

In order to run the ATLAS code, there are other device characteristics that needed 
to be defined. The first parameters that needed to be given values were the contact 
resistances. Once the resistances were defined, the bias conditions were established. 

Contact resistances were applied to the electrode contacts to simulate the 
resistances seen in the SCU experiments. The electrode contacts were placed at the very 
edge of the device as part of the gold contacts. The placement isolated the electrical bias 
from regions with potentially large temperature gradients. The resistances applied to the 
electrodes were in the low kQ range. These resistance values corresponded to the best 
electrical contacts that were achieved by the SCU group. No resistance was given to the 
CNF material but rather incorporated into the contact resistances at the electrodes. This 
is a valid assumption if the transport is considered to be ballistic [11]. 

A basic voltage sweep was determined to be the most effective bias for driving 
the device. A voltage sweep is an incrementally increasing applied voltage to an 
electrode contact over a specific range measuring a parameter. For low voltages, step 
sizes were kept the same. As the voltage increased, the different structural lengths 
reached their maximum current densities at different potentials. Therefore, the end 
voltage of the sweep was determined based on the demands of the individual structures. 
More specific information about biasing values is provided in Chapter IV under the 
relevant sections of each structure. The implementation of the bias can be viewed in the 
Appendix. 
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C. CHAPTER SUMMARY 

In this chapter, the basies of the SCU experiments and an overview of the 
eonstruetion of the simulation were eovered. From the SCU experiments, the test data 
showed that as the intereonneet length deereases, the maximum eurrent density inereases. 
Also, heat aeeumulated in the eenter of suspended regions of the deviee. 

The Silvaeo software and implementation of the SCU experimental data into the 
simulation was also eovered in this ehapter. A total of thirteen struetures were generated 
for testing. The material parameters of the CNF materials were presented and explained. 
Finally, a brief explanation was given for some of the other simulation parameters. 

In the next ehapter, the simulation data is analyzed and eompared to the SCU 
experimental data. The ehapter is broken up into three main seetions. The first two 
seetions deal with the 2D struetures’ heating eharaeteristies and eurrent densities. The 
heating in a 3D strueture is eompared to that in a 2D strueture of similar length in the last 
seetion. 
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IV. EXPERIMENTATION AND RESULTS 


Using the devices and structures discussed in Chapter III, we biased the device 
structures using a voltage sweep with the goal of determining the maximum current 
density of the simulated material. The maximum current density {Jmax) versus length with 
W-Au contacts was examined in the first group of experiments. The effect of having the 
center of the CNF in contact with substrate was examined in the second group of 
experiments. Finally, in order to verify the heating characteristics of the device, a 3D 
simulation was run and compared to its 2D equivalent. 

A. 2D SUSPENDED W-AU CONTACT 

1. Device Structures 

A total of six different lengths of CNF interconnects were created for comparison 
with the SCU experimental data. The lengths chosen included 1.25, 1.75, 2.00, 2.50, 
3.75, and 5.00 pm. These lengths were chosen to give the best range of values that 
correspond to the experimental data from SCU. The measure of length corresponds to 
the distance bridged by the CNF between the two gold contacts. It does not represent the 
physical dimensions of the CNF material. The 2D device structure was discussed in 
greater detail in Chapter III. 

2. Voltage Sweep 

Each device structure was simulated with a voltage sweep. The two 
characteristics studied included the heating and the current density within the device. 
The voltage steps were kept the same for all the different lengths of CNF. However, as 
the length increased the device reached its maximum current density at a lower voltage. 
This meant the maximum voltage was limited by the length of the device. The voltage 
was increased until a maximum current density was reached. 
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a. 


Heating of Device 


The heating in all six device structures is shown in Figures 18 through 23. 
The figures are presented in order of increasing length. In all cases, the center of the 
interconnect heated up the fastest. 



Figure 18. Ideating in a suspended 1.25 pm CNF interconnect. 
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Figure 19. Ideating in a suspended 1.75 pm CNF interconnect. 
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Figure 20. Fleating in a suspended 2.00 pm CNF interconnect. 
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Figure 21. Fleating in a suspended 2.50 pm CNF interconnect. 
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Figure 22. Heating in a suspended 3.75 pm CNF interconnect. 



Figure 23. Heating in a suspended 5.00 pm CNF interconnect 
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In the figures, the maximum temperature in Kelvins (K) varies from one figure to 
the next. This occurs as a result of different maximum voltage levels. All effort was 
made to stop the simulation near the experimental breakdown point of 900 K. 

b. Effects of Length on Jmax 

In order to measure the current density and heating throughout the 
simulation, a probe point was placed at the center of the device in the bottom conducting 
region. The data shown in Figure 24 is an illustration of the lattice temperature versus 
the current density. As the voltage increased, the current density eventually reached a 
maximum point. Each line represents a different length of interconnect. The longest 
interconnect is displayed on the far left. In general, the maximum current density 
occurred between 900 and 1000 K. This is a good sign since the device would have 
failed within this region. Any data beyond the maximum current density is unrealistic 
and should be dismissed. 


ATLAS OVERLAY 
Data from multiple files 



Figure 24. Temperature versus current density of suspended interconnects. 
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From Figure 24, it is difficult to tell where the maximum current density occurs. 
The individual curves were inspected more closely, as shown in Figure 25. In every case, 
the line had a maximum value that was represented by the bump seen in Figure 25. 
However, as just stated, any measurement above about 1000 K becomes unreliably 
because the device would have started to vaporize. The exact values obtained are 
compared to the SCU data in the next section. 


ATLAS 
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Figure 25. Close-up of temperature versus current density of a 2.00 CNF 

interconnect. 


3. Comparison with SCU Data 


The SCU data used for comparison is shown in Figure 13 in Chapter III. The 
fitted curve from this data follows [3,5], 


J 


5.39 
L ’ 


(16) 


where J is the current density and L is the length of the interconnect. In Table 6, the 
maximum current density from the fitted curve of the SCU experiments and the simulated 
data are shown. The graph shown in Figure 26 displays the data from Table 6. 
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Table 6. SCU J^ax versus simulation Jmax- 


Length (|am) 

1/L 

SCUJ^ax 

Simulation Jmax 

1.25 

0.8 

4.31 X 10^ 

3.86 X 10'’ 

1.75 

0.57 

3.08 X 10^ 

3.15 X 10^ 

2.00 

0.5 

2.70 X 10^ 

2.92 X 10^ 

2.50 

0.4 

2.16 X 10^ 

2.52 X 10'’ 

3.75 

0.27 

1.44 X 10^ 

1.93 X 10^ 

5.00 

0.2 

1.08 X 10^ 

1.56 X 10'’ 



Figure 26. Current density versus inverse length of SCU fitted data and simulation 

results. 
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When compared to the SCU data, the simulation curve has a slightly higher 
recorded current density for longer CNFs. However, the slope of the line starts to match 
up with that of the SCU data. The density of states was found to have a direct effect on 
the maximum current density without having any effect on the slope of the curve. The 
error in the results may be a result of incorrect density of states for the CNF. Next, 
supported structures are examined. 

B. 2D SUPPORTED W-AU CONTACT 

1. Device Structures 

A total of six device structures were tested for this set of data. Three different 
lengths of CNF, 1.25, 2.50, and 5.00 pm, were chosen. Each length of CNF was in 
contact with either 25% or 50% of the substrate, giving a total of six device structures for 
testing. As discussed in Chapter III, the CNF is connected to the substrate via a small 
block of SiOi material placed at the center of the CNF length. 

2. Voltage Sweep 

Just like the first set of test structures, all these devices had different maximum 
voltage values. Since the addition of the material allowed heat flow to substrate, the 
maximum voltages in the simulation needed to be increased. As in the first test, the 
voltage steps were kept the same for the simulation to create the best correlation. 

a. Heating of Device 

The heating of the six device structures are shown in Figures 27 through 
32. The majority of the heat from the devices dispersed into the substrate through the 
small block of substrate material. In most simulations the contact with the substrate in 
the center of the device caused the hottest point on the device to move towards the 
contacts. This heating effect appears more clearly for longer CNF lengths. This implies 
there is a minimum amount of unsupported material needed to create two distinct hot 
points. In Figures 29 and 30, the 25% supported device has two distinct hot points, but 
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the 50% supported device does not. In order to better understand this heating effect, two 
extra probe points were placed at the center of the unsupported regions of each device. It 
was found that, as the length of the interconnect increases, the non-centered probe points 
heated up faster but eventually reached equilibrium with the centerpoint. 



Figure 27. Heating in a 25% supported 1.25 pm CNF interconnect. 
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Figure 29. Heating in a 25% supported 2.50 pm CNF interconnect. 
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Figure 30. 


MicroTK 

Fleating in a 50% supported 2.50 pm CNF intereonnect. 
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"igure 31. 


Heating in a 25% supported 5.00 pm CNF intereonneet. 
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Figure 32. Heating in a 50% supported 5.00 pm CNF intereonneet. 
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b. Effects of Support on Jmax 

It was determined that supporting any amount of the strueture inereased 
the maximum eurrent density through the intereonneet before failure. Temperature 
versus eurrent density curves for 25% and 50% supported structures, respectively, are 
illustrated in Figure 33 and 34. The three different probe points are graphed for each of 
the structures. The shortest interconnect length is on the far right of the graph. For the 
1.25 and 2.5 |am structures, the curves completely overlap, but for the 5.0 jam structure, 
the right and left probes exceed the temperature at the center probe. A graph of the 5.0 
pim structure can be seen more clearly in Figure 35. In Figure 35, the global device 
temperature is also shown for comparison with the probes. It overlaps both the data lines 
for the left and right probes. 
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Figure 33. Temperature versus current density of 25% supported structures. 


52 




















ATLAS OVERLAY 



0 1e+06 2e+06 3e+06 4e+0€ 5e+06 6e+06 

Current_Density_Left 


Figure 34. Temperature versus current density for 50% supported structures. 
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Figure 35. Temperature versus current density for 50% supported interconnects of 

5.00 tim. 
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3. 


Comparison with SCU Data 


From the data in Figure 13, there was no equation derived for the supported cases 
of the device, which means there was no exact data for comparison. Instead, the general 
curves shown in Figure. 13 were compared to the simulation curves. The maximum 
current densities for the various structures are contained in Table 7, and the graph of the 
values is shown in Figure 36. 


Table 7. Simulated maximum current density data of supported structures. 


Length (pm) 

1/L 

25% supported Jmax 

50% supported Jmax 

1.25 

0.8 

5.06 X 10'’ 

5.714 X 10'’ 

2.50 

0.4 

3.788 X lO'’ 

4.264 X 10'’ 

5.00 

0.2 

2.847 X lO'’ 

3.446 X 10'’ 
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From just a quick look at Figure 36, the supported cases produce higher maximum 
current densities. Unlike the predicted SCU data, the curves are nearly parallel to the 
suspended case. It should be noted that the SCU group did not test any short length 
CNFs in their experiments, and as such, the curve in Figure 13 is a prediction. The 
values achieved were higher than expected but still followed the basic trend of the SCU 
data. The benefit of adding more substrate support under the interconnect appears to 
diminish as more support is added. The heating effect in a 3D structure versus a 2D 
structure is examined in the next section. 

C. 3D MODEL W-AU CONTACT 

1. Device Structure 


ATLAS 

Data from DevvoltS.str 



Materials: 

OxyNitride 
I 4H-SiC 
Tungsten 
I Gold 
□ siO-2 


Figure 37. 3D structure of a 2.00 pm CNF interconnect. 


The device structure for the 3D case differed slightly from the 2D models at the 
tungsten contact points. The contacts were placed at the very end of the CNF, which 
stretched for the entire length of the gold contacts. By placing the tungsten away from 
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the interconnect gap, the device model more accurately simulates a physical device. The 
finished structure is shown in Figures 37 and 38. The layers of CNF material and the 
contact are clearly illustrated in Figure 38. 

ATLAS 

Data from DevvoltB.str 


Y 



Figure 38. Side view of a 3D CNF interconnect. 


2. Test Parameters 

Originally, this thesis was meant to be done entirely with 3D structures; however, 
as testing progressed, it became clear that limitations in the Silvaco software would 
hinder that capability. As such, the 3D structure was biased as much as possible simply 
to show heating in the structure. All of the material parameters were kept the same and 
the device was run with a voltage sweep. 

3. Comparison with 2D model 

The resulting heating profile illustrates that the 2D model does not dissipate heat 
like the 3D model. The primary cause stems from the lack of contact between the two 
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conducting regions in the 2D model. In the 3D model, the heat ean flow around the 
eylindrieal surface of the strueture and will accumulate as far away from the substrate and 
heat si nk s as possible. The 3D profile picture is shown in Figure 39. A elose-up on the 
aetive region, in Figure 40, illustrates how the heat accumulates at the top of the center of 
the CNF deviee. 

In contrast, the 2D simulations produced a uniform heat gradient over the CNF 
materials. A close up on one edge of the hottest region of the CNF materials is illustrated 
in Figure 41. There seems to be no differenee in the top and bottom of the CNF. The top 
of the image is the side closest to the substrate of the deviee. 

The amount of material eondueting differs between the 2D and the 3D models. In 
the 3D model, about 75% of the CNF is eondueting material, whereas the 2D model has 
only 50% of the material eondueting. This would lead to more heat being generated in a 
physical device versus heat being generated in a 2D model. However, the amount of 
contact material is much greater in the 3D model, whieh allows for greater amount of 
heat to dissipate into the eontaet. It is unelear as to whieh faetor would contribute more. 



-10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 

Microns 


Figure 39. Cross-seetion view of heating in 2.00 pm 3D CNF intereonnect. 
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Figure 40. 
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Close-up of heating in a 2.00 pm 3D CNF interconnect. 


ATLAS 

Data from 2D200Lfull.str 
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Figure 41. Close-up of heating across CNF regions in a 2D simulation of a 2.00 pm 

interconnect. 
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D. CHAPTER SUMMARY 


In this chapter, the experimental data was presented and analyzed. The 2D 
simulations provided data that elosely approximated the SCU experimental data. The 
maximum eurrent density inereased as the length of intereonneets shrank. Also, the 
addition of substrate support material inereased the average maximum eurrent density of 
intereonneets. However, it was found that the heating in the 2D simulations did not 
mateh the heating within a 3D strueture. This implies there might be an error in the 
physies for the 2D model simulation. 

In the next ehapter, the eonelusions drawn from the experiments are presented. In 
addition, reeommendations for improving the eurrent model and possible future projeets 
are diseussed. 
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V. CONCLUSION AND RECOMMENDATIONS 


A. CONCLUSION 

A method to simulate the heating and maximum eurrent density for a CNF 
intereonneet was presented in this thesis. The simulation provides a first step towards 
ereating a eomplete model of a CNF intereonneet that has the ability to aeeept any input 
eharaeteristies and aeeurately model the internal physies. Further, the simulation 
demonstrated the ability to estimate the maximum eurrent density with only marginal 
error. The error present eould possibly be removed by better definition of the material 
eharaeteristies and physies of the CNF material. 

The results from the 2D simulations provided data that elosely matehed the SCU 
experiments. The most important faetor that eontributed to the maximum eurrent density 
was determined to be the density of states of the CNF material. Any error in the results 
eould be attributed to a misealeulation of the density of states. Also, the eontaet 
resistanees were found to eontribute little to the maximum eurrent density unless the 
order of magnitude ehanged. As long as the resistanees remained under 10 kQ, there was 
little to no ehange to the eurrent density eurve. 

A signifieant limitation of the model in its eurrent form is the inability to apply 
eonstant voltage or eurrent over time. If the eurrent model is driven over time, the deviee 
heats up too quiekly for the solver to eonverge. Future work would involve reetifying 
this problem and eomparing the results to the eonstant eurrent experiments done at SCU. 

As work progresses, it would be valuable to have the ability to direetly eompare 
eurrent metal intereonneets versus CNF intereonneets. As deviees shrink, Cu 
intereonneets are beeoming a greater reliability issue. With this model, it would be 
possible to study the advantages of integrating CNFs in a large seale manner. 

The next seetion eontains some reeommendations for advaneing the funetionality 
of the deviee model. Also, some issues that arose during the proeess are diseussed. 
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B. RECOMMENDATIONS 


During this project, it became quite clear that Silvaco ATLAS can be very 
limiting when creating custom materials and complex geometric structures. If one were 
to attempt a similar project that required creating a new material or complex geometry, 
ATLAS is not recommended. A wise alternative to ATLAS would have been COMSOL. 
COMSOL is a finite element solver that can simulate a wider variety of physics problems 
than ATLAS. 

ATLAS is particularly good for analyzing characteristics of well-defined 
materials, such as silicon or gallium nitride. Since there is no way to make an entirely 
new material, an existing material needs to be modified in order create a material with the 
characteristics needed. This means that all the parameters of the material need to be 
closely analyzed so that no errors occur within the simulation. 

Also, when constructing a structure with complex geometries, ATLAS takes 
considerably longer to simulate the structure. The mesh around the interfaces of complex 
geometries needs to be very well defined, which leads to longer simulation times. 
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APPENDIX. DECKBUILD INPUT CODE 


############################################################# 


# Author: Jason Brunton # 

# Naval Postgraduate School # 

# 13 September 2011 # 

# # 

# CNF interconnect Simulation file # 

# # 

# Programs: # 

# Deckbuild™ # 

# Atlas™ # 

# DevEdit™ # 


############################################################# 


go atlas 

mesh mfile=2D175L2.str conductor 

# Contact Parameters 

contact name=gate resistance=5e3 
contact name=dram resistance=5e3 

# Band Parameters 

material material=OxyNitride permittivity=5.4 
material material=4H-SiC permittivity=5.4 
material material=4H-SiC eg300=.026 
material material=4H-SiC NC300=3.0el7 NV300=3.0el7 
material material=4H-SiC NC.F=0 NV.F=0 

# Recombination Parameters 

material material=4H-SiC taun0=0 taup0=0 etrap=0 
material material material=4H-SiC augn=0 augp=0 
material material=4H-SiC augkn=0 augkp=0 
material material=4H-SiC kaugcn=0 kaugcp=0 
material material=4H-SiC kaugdn=0 kaugdp=0 
material material=4H-SiC copt=0 

# Thermal Parameters 

material material=4H-SiC tcon.const tc.const=0.12 
material=4H-SiC nsrhn=0 nsrhp=0 
material material=4H-SiC ksrhtn=0 ksrhtp=0 
material material=4H-SiC ksrhcn=0 ksrhcp=0 
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material material=4H-SiC ksrhgn=0 ksrhgp=0 
material material=4H-SiC HC.A=1.75 HC.B=0 HC.D=0 
material material=OxyNitride teon.eonst te.eonst=0.12 HC.A=1.75 

mobility material=4H-SiC munO= 100000 mup0= 100000 

models srh joule.heat fermi auger lat.temp gr.heat print 

thermeontaet num=l x.min=0.0 x.max=10.0 y.min=0 y.max=0.1 alpha=1.4e3 

method bloek newton tol.ltemp=1.0 olimit=le-4 

solve init 

solve previous 

save outlile=2D175Linit.str 

log outfde=2D 175LIVl.log 

probe name=Current_Density J.Total x=5.0 y=l.l 
probe name=Lattiee_Temp Lat.temp x=5.0 y=l.l 

solve vgate=0.001 

solve vgate=0.005 vstep=0.03 vfmal=0.095 name=gate 

solve vgate=0.10 vstep=0.40 vfinal=0.9 name=gate 

solve vgate=1.0 vstep=0.5 vfinal=7.0 name=gate 

solve vgate=7.5 vstep=0.25 vfmal=18.75 name=gate 
save outlile=2D175Lfull.str 

quit 


64 



LIST OF REFERENCES 


[1] R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical Properties of Carbon 
Nanotubes, Imperial College Press, London, 1998. 

[2] T. Saito, T. Yamada, D. Fabris, H. Kitsuki, P. Wilhite, M. Suzuki, and C. Y. 
Yang, “Improved eontaet for thermal and eleetrieal transport in earbon nano fiber 
intereonneets,” Phys. Lett., vol. 93 (10), no. 102108, 2008. 

[3] T. Yamada, T. Saito, D. Fabris, and C.Y. Yang, “Transport Phenomena in Carbon 
Nanostruetures,” IEEE International Nanoelectronics Conference, Hong Kong, 
Paper KS420, 2010. 

[4] T. Yamada, H. Yabutani, T. Saito, and C. Y. Yang, “Temperature dependenee of 
earbon nanofiber resistanee,” Nanotechnology, Vol. 21, pp. 265-707, 2010. 

[5] T. Yamada, T. Saito, D. Fabris, and C. Y. Yang, “Eleetrothermal Analysis of 
Breakdown in Carbon Nano fiber Intereonneets,” IEEE Elec. Dev. Lett., vol. 30 
(5), pp. 469-471,2009. 

[6] H. Yabutani, T. Yamada, T. Saito, and C. Y. Yang, “Joule-heating dependenee of 
earbon nanofiber resistanee,” IEEE Nanotechnology Materials and Devices 
Conference, Traverse City, MI, pp. 166-168, 2009. 

[7] H. Kitsuki, T. Yamada, D. Fabris, J. R. Jameson, P. Wilhite, M. Suzuki, and C. Y. 
Yang, “Length dependenee of eurrent-indueed breakdown of earbon nano fibers,” 
Appl. Phys. Lett., vol. 92, no. 173110, 2008. 

[8] Q. Ngo, T. Yamada, K. Gleason, M. Suzuki, H. Kitsuki, A. M. Cassell, and C. Y. 
Yang, “Eleetrothermal Transport in Carbon Nanostruetures,” Proeeedings of the 
UESI Proeess Integration 6 Symposium, Eleetroohemieal Soeiety Meeting, 
Vienna; ECS Trans., vol. 25, issue 7, pp. 487-493, 2009. 

[9] S. Maeda, T. Yamada, P. Wilhite, H. Yabutani, T. Saito, and C. Y. Yang, “W- 
deposited eontaets with earbon nanofiber using foeused ion and eleetron beams,” 
IEEE Nanotechnology Materials and Devices Conference, Monterey, CA, pp. 
125-127, 2010. 

[10] T. Yamada, T. Saito, M. Suzuki, P. Wilhite, X. Sun, N. Akhavantafti, D. Eabris, 
and C. Y. Yang, “Tunneling between earbon nanofiber and gold eleetrodes,” J. 
Appl. Phys., vol. 107, no. 044304, 2010. 

[11] T. Yamada, E. R. Madriz, and C. Y. Yang, “Induetanee in One-Dimensional 
Nanostruetures,” lEE Transactions on Electron Devices, vol. 56, p. 1834-1839, 
2009. 


65 



[12] S. lijima, “Helical Microtubules of Graphitic Carbon,” Nature, vol. 354, pp. 
56-58, 1991. 

[13] S. lijima and T. Ichihashi, “Single-shell Carbon Nanotubes of 1-nm diameter,” 
Nature, vol. 363, pp. 603-605, 1993. 

[14] M. S. Dresselhaus, G. Dresselhaus, and R. Saito, “Physics of Carbon Nanotubes,” 
Carbon, vol. 33, no. 7, pp. 883-891, 1995. 

[15] Q. Ngo, T. Yamada, M. Suzuki, Y. Ominami, A. M. Cassell, J. Li, M. 

Meyyappan, and C. Y. Yang, “Structural and Electrical Characterization of 
Carbon Nano fibers for Interconnect Via Applications,” IEEE Trans. Nano., vol. 6, 
no. 6, pp. 688-695, 2007. 

[16] C. V. Thompson, “Carbon nanotubes as Interconnects: Emerging Technology and 
Potential Reliability Issues,” Reliability Phys. Sym., pp. 368-376, 2008. 

[17] J. Ei, A. Cassell, H. T. Ng, R. Stevens, J. Han, and M. Meyyappan, “Bottom-up 
approach for carbon nanotube interconnects,” Appl. Phys. Lett., vol. 82, no. 15, 
pp.2491-2493, 2003. 

[18] N. Srivastava, R. V. Joshi, and K. Banerjee, “Carbon Nanotube Interconnects: 
Implications for Performance, Power Dissipation and Thermal Management,” 
Electron Devices Meeting. lEDM Technical Digest, pp. 249-252, 2005. 

[19] N. Kuruvilla and J. P. Raina, “Statistical Eatency Analysis of Carbon nanotube 
Interconnects due to Contact Resistance Variations,” Inter. Conf. 
Microelectronics, pp. 296-299, 2008. 

[20] A. Srivastava, Y. Xu, and A. K. Sharma, “Carbon nanotubes for next generation 
very large scale integration interconnects,”/. Nanophotonics, vol. 4, no. 041690, 
2010 . 

[21] Silvaco, Inc. Atlas user’s manual. Santa Clara, CA: Silvaco, Inc., 2010. 


66 



INITIAL DISTRIBUTION LIST 


1. Defense Teehnieal Information Center 
Ft. Belvoir, Virginia 

2. Dudley Knox Library 
Naval Postgraduate Sehool 
Monterey, California 

3. Engineering and Teehnology Currieula Offiee, Code 34 
Naval Postgraduate Sehool 

Monterey, California 

4. Chairman, Code EC (590) 

Department of Eleetrieal and Computer Engineering 
Naval Postgraduate Sehool 
Monterey, California 

5. Professor Todd R. Weatherford, Code EC (590)/Wt 
Department of Eleetrieal and Computer Engineering 
Naval Postgraduate Sehool 

Monterey, California 

6. Professor Sebastion Osswald, Code PH 
Naval Postgraduate Sehool 
Monterey, California 

7. Jason A. Brunton, Code EC (590) 

Department of Eleetrieal and Computer Engineering 
Naval Postgraduate Sehool 
Monterey, California 


67 



